Loss of Andreev Backscattering in Superconducting Quantum Point Contacts 
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We study effects of magnetic field on the energy spectrum in a superconducting quantum point 
contact. The supercurrent induced by the magnetic field leads to intermode transitions between the 
electron waves that pass and do not pass through the constriction. The latter experience normal 
refiections which couple the states with opposite momenta inside the quantum channel and create 
a minigap in the energy spectrum that depends on the magnetic field. 
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The transport in superconducting-normal-mctal (SN) 
mesoscopic hybrid structures is controlled by normal and 
Andreev reflections. The fundamental property of the 
Andreev reflection is its almost exact backscattering: tra- 
jectories of an incoming particle and reflected hole coin- 
cide within an angle of a maximum {^fO"^ where 
^ is the superconducting coherence length and kp is the 
Fermi wave vector. Backscattering holds also in the pres- 
ence of magnetic field since the electron-hole conver- 
sion length ^ ^ is much smaller than the Larmor ra- 
dius tl, i.e., £,/rL ^ {B / Hc2){kF£)~^ ■ This has been 
checked, for example, by the electron focusing technique 
0. Deviation from backscattering come from interfer- 
ence of electron and hole waves with an inhomogeneous 
order parameter phase. Effects of a supercurrent-induced 
transverse force on backscattering have been considered 
in Refs. 0,0. In Refs. the phase difference introduced 
between the partial waves resulted in a decrease in the 
reflected hole intensity leading to an interplay between 
the normal and Andreev processes. 

Though deviation from exact backscattering is small, 
it can still be noticeable if it is comparable with the size 
of the setup. A convenient device that satisfies this re- 
quirement is a ballistic superconducting quantum point 
contact. In the present paper we study how the vi- 
olation of the fundamental semiclassical backscattering 
property of the Andreev reflection affects the subgap en- 
ergy spectrum in such contact. Under the condition of 
exact backscattering, each electron that passes through 
the constriction is reflected as a hole that returns along 
the same trajectory and thus has also to pass through 
the constriction. The energy spectrum of subgap states 
is then E = iApCOSx where Aq is the superconduct- 
ing gap and 2% is the phase difference between the elec- 
trodes 0, Q ■ The two energy branches which correspond 
to quasiparticles propagating through the constriction in 
opposite directions cross the Fermi level = at 2x = tt. 
We show that the loss of exact backscattering leads to a 
dramatic change in the spectrum such that a minigap 
appears near 2% = tt. The deviation from backscatter- 



ing can be produced, for example, by an exchange, dur- 
ing the Andreev process, of a Cooper pair momentum 
induced by an applied magnetic field. This momentum 
mixes the channel modes with the modes that do not pen- 
etrate inside but are normally reflected from the channel 
end. The normal reflections couple the waves propagat- 
ing through the constriction in the opposite directions 
and lead to formation of a minigap in the energy spec- 
trum similar to that for contacts with normal scatterers 
Varying the magnetic field one can tune the 
degree of normal reflection and manipulate the minigap 
thus controlling the transport properties of the contact. 

Model.- The loss of exact backscattering at the An- 
dreev reflections in a quantum point contact can be more 
clearly illustrated for an auxiliary structure shown in Fig. 
n (a): A single mode channel with a radius a ^ kp^ is 
open into a normal semi-spherical region with a radius R 
much larger than the superconducting coherence length 
^. The normal region is surrounded by a superconduc- 
tor which carries a supercurrent with a momentum ftk^ 
perpendicular to the channel axis. For R ^ ^ one can 
describe quasiparticle propagation using a trajectory rep- 
resentation. Due to the transfer of Kkg the Andreev re- 
flected trajectory deviates from its initial direction 
such that it can miss the constriction and experience nor- 
mal reflections from the insulating barrier. The trajec- 
tory returns to the constriction after several reflections, 
which results in coupling of states propagating through 
the constriction in the opposite directions. The transfer 
of ks causes a trajectory deflection by an angle kg/kp, 
thus the shift of the trajectory over a distance R would 
be kgR/kp. The probability of normal reflection thus 
depends on the ratio of the trajectory shift to the trans- 
verse channel dimension a. For a single- mode quantum 
channel a ~ kp^, this ratio is kgR/kpa ^ k^R. 

In a superconducting point contact, the quasiparti- 
cle wave functions for subgap states decay over dis- 
tances of the order of i? ^ ^, which gives an estimate 
ksR ~ kgS, < 1. However, the trajectory shift ks£,/kp 
is less than the wave length so that the trajectory pic- 
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FIG. 1: (a) A single-mode channel is open to a normal region 
(white semicircle) in a contact with a superconductor (grey 
region). Andreev reflected trajectories deviate from initial 
direction due to the transverse pair momentum fika and ex- 
perience normal reflections from the insulator surface (black) , 
which couple right-moving and left-moving u~ states, (b) 
Asymmetric and (c) symmetric point contacts. 



ture is not applicable. Moreover, a single-mode contact 
with a radius a ^ kp^ radiates an electronic wave which 
is determined by diffraction. To match the particle 
and hole waves propagating along the channel with 
the quasiparticle waves in the corresponding supercon- 
ductors we assume, for simplicity, that the right super- 
conductor occupies a halfspace z > and has a specularly 
reflecting surface. We introduce the spherical coordinates 
X = r sin 6 cos 4>,y — r sin 9 sin (p, z = r cos 6 and expand 
the wave functions in spherical functions with odd angu- 
lar momenta. This ensures vanishing of the microscopic 
wave function at the superconducting/insulator bound- 
ary 6 = tt/2. Far from the origin (r ^ a) we have 



(1) 



Different angular harmonics are orthogonal to each other 
within < 9 < n/2. For a waveguide a < kp^, the 
radiated/incident diffraction field is expi^iikpr) cos9/r. 
We now assume that it is this only mode proportional 
to Pi.o = cos 9 that ideally transforms into the channel 
mode Mo, vq without reflections, while all other modes I ^ 
1 are normally reflected from the waveguide end without 
transmission into the channel. 

Scattering Matrix.- For more quantitative description 
we use the Bogolubov - de Gennes (BdG) equations: 



, ihV--AY-EF 
2m V c / 



, ihV + -aY - Ep 

2m \ c / 



u + Av = Eu , (2) 

V - A*u = -Ev , (3) 



where A is the vector potential of the magnetic field 
B = B{z)x, and A is the gap function. We assume a 
step-like form A = Aoe*'-'^+''='"\ where xr.l = ±X is the 
zero-field constant phase in the right (left) superconduc- 
tor with 2x phase difference between the superconduc- 
tors; ks is a constant wave vector. It enters the supercon- 
ducting velocity mVg = ^ik^ — (2e/c)A that determines 



the difference in the eikonals of particle and hole wave 
functions, (m/h) J Vg - dr. The magnetic field is screened 
in the superconductor, A = ALi3oy(exp(— z/A^) — 1), 
at the London length A^. From the screening condition 
Vg = at large z we find = (27r/^o)ALi?oyj where (/)o 
is the flux quantum. Assuming for simplicity Xl ^ £, we 
neglect A in the region r ~ ^ where the wave functions 
are localized. The parameter kg^ that determines the 
probability of normal reflections of the channel modes is 
ks^ ~ Bo/Hcm where Hem ~ 0o/(Al^) is the thermo- 
dynamic critical fleld of the superconductor. The gap 
Ao is suppressed near the superconductor surface. How- 
ever, this does not change the backscattering properties 
of Andreev reflection; we ignore it in what follows. 

Inside the single-mode channel there are two particle 
and two hole waves oc e^^^^^ with amplitudes Uq and 
v^, respectively, corresponding to the momentum pro- 
jections ±hkz on the z axis. A particle Uq and a hole Vq 
propagate in the +z direction while a particle Uq and a 
hole Vq propagate in — z direction. Using the scheme em- 
ployed in Ref. we introduce the scattering matrices 
Sni^, Xi ^sr) and Sl{£, —X, ^sl) that relate the incident 
and outgoing wave amplitudes respectively at the right, 
z = d/2, and left, z — ~d/2 ends of the channel: 



= Sr 



(4) 



Here d <C ^ is the channel length, ks/, and Ilsr are the 
superflow momenta in the left and right superconductors, 
respectively. The wave functions at the both ends of the 
channel have different phase factors: 



-"0 / R \ "0 / L 

The solvability condition of Eqs. Q and |5j yields 



(5) 



det 1 



.^^^k^dg^^za^k^dg^ ) = . (6) 



The matrix S is unitary: SS'' = 1. Indeed, the BdG 
equations (01, ® conserve the quasiparticle flow 



div 



(^-ihV - -Aj u + u (^ihV - -A ) M* 



ihV 



V — V (^i 
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Since this flow vanishes for \E\ < Aq deep in the super- 
conductor it should be zero also in the channel, whence 
Wol^ + kcTP = 1^0 P + \vo\^ which results in SS'< = 1. 

We now calculate the matrix S explicitly with the ac- 
count of reflections from the channel end. Wave functions 
decaying into the right superconductor obey the relations 



'a+(ks_R)e 
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(7) 
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that couple the electron and hole amplitudes near the 
channel end |r| ^ ^. Here 



^ — hvp/ ^0, and e = E/Aq. For the left superconductor, 
similar expressions hold with x ~* ~X and a"*" — a~. 

Let us consider the matrix S at the right end of the 
channel. We place the origin of the coordinate system at 
the right channel end and consider the case x = since, 



generally, S{x) 



S(x = 0)e 



According 



to Eq. J^l, the wave function amplitudes in the right 
superconductor at distances a <C r <C ^ from the origin 
have the form (we omit the index R) 

U+ = u+P,,o + , V+ = v+P,,o + K , 
U- = Uq Pi^o + , V- = VPi,o + 'i'v ■ 

The amplitudes '^^ and '^y stand for the sum of com- 
ponents with ^ 7^ 1 in Eq. Q and describe the modes 
which experience normal reflections at the channel end. 
Introducing the operator that describes this reflec- 
tion we can write the relation between the amplitudes 
^+ = Re'^u, = R-e'^v The functions ^u, *u are or- 
thogonal to the mode Pi^: (Pi,o|^u,'u) = 0. The angular 
brackets denote the angular average within < 9 < tt/2. 

The modes and ^„ together with the traversing 
mode mq, vq experience Andreev reflections while only 
Mo and vo contribute to the flow through the channel. 
The unitarity SS^ = 1 implies that the quasiparticles 
which are scattered normally off the superconductor sur- 
face and the channel end will eventually return into the 
constriction either as particles or as holes after certain 
number of Andreev reflections at the superconducting 
side. The Andreev relations (|7I8() should be applied to 
the total functions U and ^ at a hemisphere with the 
center at the origin and the radius much smaller than ^ 
where kjr ~ fc^rsin^sin^ <^ 1. Taking the derivatives 
only of the rapidly varying exponents in Eq. we obtain 



= e-"^^(M+Pi,o + i?e*«) , (9) 
- -e''^-=(VPi,o + ^'.) , (10) 



where 



= e - ^k,r/2 + i^l - (e - efcsr/2)2 



(11) 



and ksr — kg sin 6 sin (f>. Since the normal reflection at the 
channel end is associated with the momentum transfer of 
the order of hkp one can neglect the dependence of R on 
energy on the scale Aq and assume R = —e^f'' where 
tfr is a constant phase shift. We then solve Eqs. 
(|10|l for the functions ^'k,^'^ and use the orthogonality 
{Pi.o\^u,v) = 0. This yields two equations that couple 
the amplitudes Uq , Vq to Uq , v'^ through the matrix 



S 



|c,|2-l)-(c,-c:)e'^'^*^a,) (12) 



where 



(Pi,o (1 + e 



Pi 



{Pi,o {e"fi<^ + e'f-') 'Pi^o) 



(13) 



One sees from Eq. (HTJ that e*'^'(fcs) = -e"^"^-^ (-fc^) 
whence c*{ks) — — c_g(— fc^). Moreover, Ce(fcs) is an even 
function of ks because a change in the sign of ks can be 
compensated by the shift (f> tt + (j> in the integral over 
the angles. Thus, c*{ks) = ~c^^{ks). For small ks£,, 



i{ksiY{PlQ sin^ ^sin^ 0)e2''' 
8(p2o) sinS 



(14) 



Without kg one has = e"' where e"' — e + iy/l — e^. 
As a result, the diagonal components of S vanish, thus 
the +pz and —pz states are decoupled. 

In the diffraction picture, the transitions that couple 
the penetrating and non-penetrating modes are caused 
by the angle-dependent Doppler shift of energy propor- 
tional to ksr in Eq. Ijllll . As a result, the wave fronts of 
reflected holes are distorted as compared to those of in- 
cident particles. The interference of these waves near the 
channel end results in the suppression of the amplitude 
of the Andreev reflected wave entering the channel. 

Results - Consider first the zero-bias conductance of a 
normal-metal - quantum-channel-superconductor junc- 
tion [13 G. = (eV^^)(l - l-^iil' + 1^121') where \Sn\^ 
and 1 5*121^ are probabilities of normal and Andreev re- 
flection, respectively. We get for small fcs^ 
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Here we introduce a fleld He ^ Hem through 
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I sm 
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(15) 



Consider now an asymmetric structure that consists 
of a superconducting tip with a curvature radius smaller 
than Ai in a contact with a bulk superconductor, see Fig. 
n(b). In this case k^^ = while k^^ = k^ ^ 0. On the 
right end of the channel the matrix Sn = S{e,Xj'^s) is 
determined by Eq. ((T^ . On the left end the matrix Sl = 
S{e,—x,0) assumes an Andreev form 5'l = e~^^e~^^'^' ax- 
The phase shift kzd — cpr drops out and Eq. © yields 
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(1 
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2(c: 



I cos(2x) . (16) 



For fc^ = and — e^^ we obtain a standard gapless 
expression e = ±cosx- For a nonzero kg, a gap opens 
in the energy spectrum. Indeed, consider Eq. H16|l in the 
limit of small magnetic fields and energies. It becomes 



1 



= cos^ X + o 



{ic, + 1)' + [icl - if 



e=0 
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Within the leading terms in B/Hc we find 



(17) 



where the minigap in the spectrum is = j{Bq/ Hc)'^ . 

In the case of a symmetric contact shown in Fig. 
n (c) the solution of the screening problem yields k^^ — 
— ksfl = — ks. The spectral equation jnj with Sr — 
S{e,Xi'^s) and Sl = 'S'(e, — x, — ks) reduces to 



cos^ X 



siv?{kzd — ipr) ■ (18) 



The spectrum has a gap in the presence of kg . In the limit 
of low energies and small magnetic fields the right hand 
side of Eq. ((TH|l can be treated as a perturbation. We put 
(|c,|2 - 1)2 ~ {Bq/HcY where is defined in Eq. HSJ 
while (ce-c*)2 « -4. At the same time, (ce+c*)2/4 ~ . 
Finally we get Eq. (|17|l where 



in the region eV < eg [lj|. Varying the magnetic field 
one can thus observe a transition from ballistic to high- 
resistance behavior of the contact. 
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Discussion.- Since the wave vector kg ^ is much 
smaller than kp it induces transitions only between the 
modes with close transverse quantum numbers. Thus, 
the predicted effect can be more easily seen for a con- 
tact that is transparent only for a few modes. On the 
contrary, for a multi-mode channel, the coupling to the 
refiected modes that mixes p and — p states has a small 
weight while the transitions occur mostly between the 
penetrating modes. For a large area SNS junction these 
transitions result in the subgap spectrum instability with 
formation of energy bands [jj]. 

Equation H17(l coincides with the spectrum in the pres- 
ence of normal scatterers fio", 11]. Note that the gap in 
a symmetric contact Eq. (|19|l vanishes for certain phase 
difference k^d — (pr — 7m. This is a result of resonant 
tunneling through a system of two barriers with equal 
reflection coefficients. The transmission probability jTp 
through such system is unity at the resonance such that 
the gap eg = 1 — |Tp disappears. A small asymmetry 
in the scattering conditions removes the resonant tun- 
neling effects so that a gap will exist for any phase shift 
kzd — (fir- The extreme asymmetric case is illustrated by 
Eq. H17(l . Similar effects of resonant tunneling and mini- 
gap oscillations as functions of k^d can also take place for 
other mechanisms of normal reflection such as interface 
barriers, mismatch in the material parameters, or small 
normal scattering from the step-like gap potential |0| . 

The predicted spectrum can be tuned by varying the 
magnetic field. The minigap is not small and can reach 
values of the order of Aq for Bq ^ Hem- It can be mon- 



itored by measuring the Josephson critical current that 
decreases with the minigap [llj. Moreover, the minigap 
affects the dynamic properties of the point contact. In 
particular, it is responsible for suppression of the time- 
averaged quasiparticle current for voltage biased contacts 
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